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Abstract 

The formalism for modeling multiple fermion generations in a warped extra di- 
mension with a soft-wall is presented. A bulk Higgs condensate is responsible for 
generating mass for the zero-mode fermions but leads to additional complexity from 
large mixing between different flavors. We extend existing single-generation analy- 
ses by considering new special cases in which analytical solutions can be derived and 
discuss flavor constraints. The general three-generation case is then treated using a 
simple numerical routine. Assuming anarchic 5D parameters we find a fermion mass 
spectrum resembling the standard model quarks and leptons with highly degenerate 
couplings to Kaluza-Klein gauge bosons. This confirms that the soft-wall model has 
similar attractive features as that found in hard-wall models, providing a framework 
to generalize existing phenomenological analyses. 
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1 Introduction 



The warped extra dimension provides an alternative framework in which to address the 
Standard Model (SM) gauge hierarchy problem pQ and the fermion mass hierarchy [2j [3] . It 
is a particularly attractive scenario because, by the AdS/CFT correspondence |4], the five- 
dimensional (5D) framework is dual to a four-dimensional (4D) strongly-coupled conformal 
field theory. This allows the physics of the warped fifth dimension to be recast in terms of 
4D strong dynamics. In particular, an infrared (IR) brane used to generate the Higgs cutoff 
scale is interpreted in the 4D dual as the breaking of conformal symmetry at low energy 
and the generation of a mass-gap by an operator of infinite scaling dimension. While the 
essential physics is captured in this theoretically idealized situation, it is more natural to 
expect operators of finite dimension in the dual theory. This can be achieved in soft-wall 
models in which the "hard- wall" IR brane is replaced by a scalar field (the "dilaton" ) whose 
nontrivial bulk profile corresponds to conformal symmetry breaking in the dual 4D theory. 

The Standard Model in the soft- wall warped dimension was considered in Ref.[5]. Since 
there is no IR brane, SM fields are necessarily bulk fields, which includes not only the gauge 
bosons and fermions, but also the Higgs field. Even though the fifth-dimension is semi- 
infinite, the dilaton does provide a dynamical cutoff to the warped dimension. This leads to 
a discrete Kaluza-Klein (KK) mass spectrum, but with the distinctive feature that there is 
a variety of KK spacing between the resonances, including linear Regge-like behavior as in 
QCD. A bulk Higgs condensate is responsible for breaking electroweak symmetry but causes 
the analysis of fermions to be particularly involved. An analytical solution can nevertheless 
be obtained in the case of a single fermion generation confirming that the nice features of 
hard-wall scenarios, such as fermion mass hierarchies and universal KK gauge couplings, also 
exist with the soft wall [3]. However these features have yet to be confirmed in a complete 
three-generation soft-wall model. 

In this paper we present a more comprehensive analysis of bulk fermions in a soft-wall 
warped dimension (see also [3, [7]). While new analytical solutions are found for special 
cases in the case of a single generation, the complete three-generation case can only be 
treated numerically. This is because the bulk Higgs condensate causes large mixing between 
fermion flavors in the equations of motion which makes finding analytical solutions nontrivial. 
Nonetheless, numerical techniques can be used and we present a numerical routine that can 
solve the general problem with arbitrary 5D mass parameters. Importantly we find that 
starting with "anarchic" 5D parameters we are able to generate fermion mass hierarchies 
and universal couplings to KK gauge bosons, analogous to that found in hard-wall models. 
This provides a framework to generalize existing phenomenological analyses to include the 
soft-wall Standard Model. 

The layout of this paper is as follows. In Section [21 we review the setup needed to 
model fermions in the soft-wall background and present the fermion equations of motion. 
In Section [31 we develop the tools needed to solve the equations of motion in the special 
cases where the equations can be partially decoupled. One of these cases has been detailed 
before in Refs. 016], while the remaining cases are new. We then show that these additional 
analytical cases can be used to reproduce many of the recent numerical results of [7\. We 
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conclude this section by discussing flavor changing neutral current processes and by detailing 
the couplings of fermions to gauge bosons-an analysis which is easily generalized to other 
bulk couplings. In Section HJ we present a very simple, non-iterative routine which can be 
used to analyze multiple generations of fermions in an arbitrary background. We present 
the full dependence of SM fermion masses on the 5D bulk mass parameters and compare 
the results to a typical hard-wall model. The behavior is shown to be very different in 
the phenomenologically interesting region of the parameter space, where the bulk SU(2)l 
doublet and singlet fermions have opposite bulk masses. We then present results for the 
case of three fermion generations with substantial mixing between bulk profiles, and find 
example spectra resembling the up- and down-type quarks (and charged leptons) in the 
spirit of Ref.[8]. 

2 Fermions in the Soft- Wall Background 

We work in a 5D spacetime z) with conformal coordinate z and metric: 

] llMN<- 



ds 2 = e~ 2A ^VMNdx M dx N 



where rj MN = diag(— ,+,+,+, +). In particular we will consider a pure AdS metric, i.e. 
A(z) = log kz with k the AdS curvature scale. The spacetime is defined on the interval 
z G [,2o? oo), where zq is the location of the ultraviolet (UV) brane. Though the spacetime 
extends to z — > oo, we have in mind a soft- wall setup in which the dilaton, $ obtains a 
background value and provides a dynamical cutoff to spacetime along the fifth coordinate z. 
In this scenario, gauge and matter fields are described by the action, 

-* 



S = j <Pxy/=ge-*£, (2) 

where £ is the 5D Lagrangian. While much of our discussion of fermions is valid in general, 
for the sake of concreteness we will specifically consider a dilaton profile given by: 

= M 2 , (3) 

where the soft-wall mass scale \x ~ 1 TeV. This form for the dilaton was also studied in 
detail in Ref.[5]. 

Consider 5D Dirac fermions, {^r) which transform as a doublet (singlet) under 
SU (2) l- It is straightforward to embed our setup in a theory with a bulk custodial SU (2) l x 
SU(2)r symmetry, but this will not be essential for our discussion. In the absence of Yukawa 
interactions, the fermion action is given by: 



S = — d xy/^g e 



(4) 
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where = e A 8^ is the vielbein and Dm = 9m + ^>m is the covariant derivative with spin 
connection uom- The index a is an SU{2) label, while i,j are 5D flavor indices. 

The projections of the Dirac spinors are given by \P£± = ±7 5 \l/2± an d similarly for ty l R . 
Dirichlet conditions are imposed on the fields and at the UV boundary: 



V R+ (x,z) 



0. 
0. 



(5) 



Without bulk Yukawa interactions these boundary conditions give rise to massless chiral 
fermions from the 4D point of view. These zero-modes can obtain a mass by introducing a 
Yukawa coupling to the Higgs, whose vacuum expectation value (VEV) is z-dependent. The 
Yukawa interaction contribution to the action is: 



5- 



A*- 7 ' 

) nl.a.ru = ~ \ (l ■'' V (x, z) H" (x, z) V R (x, Z) + h.C. 



Sx^ge'^ m ij (z) ^ L {x, z)^ R (x, z) + h.c. 



where we have substituted the background value for the Higgs field: 

h(z) 



H(x,z) 



Hz) 



V2 



(6) 



(7) 



and dropped the SU(2) labels *f?L = ^ 2 L . The effective ^-dependent bulk mass term arising 
from the Yukawa interaction is simply: 



m 



z = 



J2k 



h(z). 



To ensure a discrete spectrum of fermion masses, the Higgs VEV must grow faster than the 
metric factor, e A ^ = l/(kz), decays. Namely]]] 



lim 



h(z) 



oo. 



z— >oo z 

Varying the action with respect to r, we find the equations of motion: 



x m lj ij) RT 



0. 
0. 



(9) 



(10) 
(11) 



where we have defined \I/ = e 2A+ ®/ 2 ip. This transformation shows that the fermion mass 
spectra do not depend on the presence of the dilaton. Rather, it is the Higgs VEV that sets 
the fermion spacing, in contrast to the case of bosonic fields. 



1 Other possiblities may also be considered. For example, if lim h(z)/z — > [i > 0, there can exist discrete 

z — >oc 

low- lying modes with a continuous spectrum above a "mass gap", as in Refs. [SJE]. 
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The KK expansion for the fields ipL,n± is assumed to be: 



^i ± {x,z) = ^/2 n) W^ (n) (x), (12) 

n,a 

^ R± (x } z) = £/g B) (z)i£ (B) (s), (13) 

n,a 

where 'y^d^ip'^^ = — m"?/>^ ri ' ) (no sum over a). Similar to the conventions of Ref.|Sj we have 
introduced separate Latin and Greek indices labelling the 5D and 4D flavor, respectively. 
Defining the vectors: 

ria(n) 



f± W = ^ J f f(n) J , (14) 

allows the equations of motion for the 5D fields to be written in the form: 

[±d 5 5« + M ij ] f ± a{n \z) = m%f^ n \ (15) 

where the mixing matrix is defined as 

-a f Ml j m ij (z) \ 
M = e A ( l,, J,' . (16) 



Note that a, i,j run from 1, . . . Np, where Np is the number of fermion generations. Thus 
M. l i is a x 2Np matrix, and Eq. (|15|) represents a coupled system of ANp differential 
equations for each a. The 4D fermion fields il)^ n \x) are canonically normalized by requiring 
that: 



dz {ftt^ fi± m) + (/S n) ) f /S m) =S nm S ap . (17) 

Note that the index i is to be summed over in this expression. This completes the general 
discussion of the fermion setup. To obtain the spectrum of fermion masses, the equations of 
motion (|15p are solved subject to the boundary conditions (jSJ) and orthonormality conditions 

m 



3 Fermion Spectrum 

The coupled equations (|T5l) cannot be solved analytically except for a few special cases, 
depending upon the particular form of the Higgs VEV and the relative bulk masses for the 
fermions. Up to this point, the only solvable cases that have been presented in the literature 
have involved just a single generation with degenerate bulk masses for the fields ^ l and ^ r. 
As it turns out, in AdS space there are additional special cases which allow for the second- 
order equations of motion to be diagonalized and solved exactly. The collection of solvable 
cases provides a qualitatively complete picture of fermion behavior in the entire parameter 
space. 
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Next, we review the single generation case in detail. We begin very generally, emphasizing 
that these methods apply to a wide variety of soft-wall models in AdS. We then specialize to a 
quadratic VEV and solve the equations of motion directly for the special cases. The analytic 
solutions allow us to verify the results of the numerical routine we present in Section H] (as 
well as a recent numerical treatment in which the Yukawas are treated perturbatively [7]). 



3.1 Single Generation 

For a single generation of fermions, equation ffl5l) becomes: 

(±d z + M)f£\z)=m n f^\z), (18) 
where M. is a 2 x 2 mixing matrix: 

M = e ' A ( M i\ ™m ) < 19 > 

^ m(z) M R J 

The equations for and /_ can be decoupled by deriving a second-order equation from 
(fTSj) . The fields /± obey a Schrodinger-like equation: 

(-d 2 z +V ± )fi n) =mlfi n \ (20) 
where the "potentials" are given by: 

V±(z) =M 2 tM'. (21) 
The difficulty in solving (|2"U1) is due to the fact that the mixing matrix generally cannot 

(n) 

be diagonalized through global transformations of the functions, fl R ±- However, there are 
special cases for which the second-order equations can be decoupled further. They occur 
whenever: 

Ml = Mr, "degenerate" 
M L + M R ± d z e A{z) = 0. "split" (22) 

The "degenerate" case is separable in any background. The "split" cases are separable 
regardless of the Higgs VEV in AdS, where the split-case condition simply becomes Ml + 
M R ± k = 0. 

For generic forms of the Higgs VEV, it is most useful to work with transformed fields, 



In this basis the equations of motion are given by, 

(±d z +M)gt\z)=mJ™\z), (24) 



where 



fi_ u t MU _£±( 2m{z)+M L + M R M L -M R , 

M-UMU- 2 I Ml -M r - 2m (z) + M L + M R »' U " 



while the boundary conditions (jSJ) become: 

We may also define transformed potentials, V±, in direct analogy with (j2l|) . For the degen- 
erate case, both of the potentials V + and V_ are simultaneously diagonal in this basis. In the 
split cases, only one of the potentials V± will be diagonal. After solving the corresponding 
pair of decoupled second-order equations, the first-order equations (1241) can then be used to 
generate the remaining solutions. 

Below, we consider the degenerate case and one of the two split cases, Ml + Mr + k = 0, 
assuming the following form for the Higgs VEV: 

h(z) = r/fc 3 / 2 /iV, (27) 

giving m(z) = bk(fxz) 2 where b = A5^/v^2, as in [5]. We also parameterize the bulk masses 
in units of the AdS curvature, Ml,r = CL,Rk, where cl,r are dimensionless coefficients. 

3.1.1 Degenerate Bulk Masses 

The solution to the degenerate bulk mass case cl = cr = c was presented in detail in [5]. 
The lowest-lying mode was found to be: 

f /, 2 2\-l/2+|c| for U > 1 

2 _ ] r(-l/2+|c|) Z 0) tOT l C l > 2 ' (2H) 



c < 



TrseccTT— 0(l/2-c)— 0(l/2+c) w l°l ^ 2 ' 

where ip is the digamma function. Note that the lowest-lying mode is very light for \iz§ <C 1 
and |c| > 1/2 (becoming exponentially small with increasing |c|), while for |c| < 1/2 the 
fermion mass is of order b\j? . 

3.1.2 Split Localizations 

Here we consider one of the "split" cases, cl + cr + 1 = 0. The other case, cl + cr — 1 = 0, 
is very similar. With this choice, the transformation (1231) will diagonalize the potential V + 
for any Higgs VEV in AdS. However, for our choice, h(z) ~ z 2 , the untransformed potential 
V- happens to be diagonal, 

c(c-l) i A2,,4^2 



b 2 ^Z 2 



so we will work in this basis. 

A consistent solution requires that either f { L n } = or f R n } = 0. This is a peculiarity of 
the particular choice of the Higgs VEV and will not be true for other forms. The result is 
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that the full tower of orthogonal solutions is most easily described in terms of two "distinct" 
KK towers of solutions. The first solution is: 



(30) 



= (32) 

/£?(*) = 0, (33) 

where iv£_ is a normalization constant and U(a, b, y) is the Tricomi confluent hypergeometric 

function. For this tower, the boundary conditions /#+| zo = and f^l\ — are equivalent. 
There is only a single orthonormality condition, 

dzff}ft ] =S nm , (34) 

which in fact implies the correct orthonormality condition for the remaining fields, 



An) Am) , An) Am) 
J L+J L+ J R+J R+ 



(35) 



(36) 



dz 

The other KK tower is given by: 

ftM = (37) 

/2W = -^(^/S (38) 

/£!(*) = 0. (39) 

where Nj^_ is a normalization constant. The spectrum for this tower is found by imposing 
the boundary condition | = 0, while the normalization condition is, 



dzf^ff =5 nm . (40) 

The lowest lying mode of the second tower is very light and approximating the mass re- 
quires some care. For wig <C 6/i 2 , we can expand the functions using techniques of so-called 
boundary perturbation theory of quantum mechanics [TOl [TT] . The function obeys the 
Schrodinger-like equation, 



(c + l)(c + 2) 2 



2' 



= (41) 
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For small tuq, we can write as a product of the zero- mode solution and a correction: 



= ((z)F(z), 



where ((z) satisfies the zero-mode equation 

-dlC 4 

The solution may be written as, 



C+1) ( C + 2) +6W 



C = o. 



«z)=N®z L '*K v W*'/2), 



(42) 



(43) 



(44) 



where iV^°2 is a constant, K v is the modified Bessel function and v = 3/4 + c/2. The function 
F(z) obeys the second-order equation, 



[C- 2 d z (( 2 d z )+ml] F(z) = 0, 



and may be expanded in powers of ?71q as, 



fSl * a*) 



pz POO 

+ m 2 dz'C 2 / dz"( 2 + 0(r4) 

J ZQ J z' 



(45) 



(46) 



Such an expansion has also been used in Ref. [12J to approximate wavefunctions in soft-wall 
models. In contrast, here we are using it to solve the boundary value problem. The UV 
boundary condition, 

' f { R ' ' 



(0) 
R+ 



■AO) 
~Jr- 



0. 



results in the following approximate expression for m^: 



m ~ 



K u (b^z 2 /2)l(z ) 



(47) 



(48) 



where 



poo 

l(z) = / dz'Ciz'). (49) 

J z 

The expression (1481) can now be expanded for small zq. For c > — 1/2 we find, 

'b^z^ c - 1/2 



2 _ 2(c/2 + 3/4) (b 2 ^z 2 



T (c/2 + 3/4) V 4 



T (-c/2 + 1/4) + r (c/2 - 1/4) 



(50) 



In the limit c ^> 1/2, this expression simplifies further to 



~2 




This expression reveals a lower bound on the fermion mass in this region of the bulk mass 
parameter space, m min = (fiz )bfi. 

For c — 1/2, the above expansions are poor approximations because the mass becomes 
0(bfi). To deal with this regime, we can apply mathematical techniques from "supersym- 
metric quantum mechanics" to determine the mass [13J. Consider a quantum mechanical 
system for which the Hamiltonian may be factorized as: 

[-d z + W{z)\ [d z + W(z)\ $ = m 2 J,. (52) 

The "superpotential," 

W(z) = + b[i 2 z, (53) 

gives rise to the "ordinary" potential for the function ip: 

Viz) = W 2 -W'= (c+1) l C + 2) + b 2 ^z 2 + (2c + l)bfi 2 . (54) 

z 2 

It is clear from (152]) that there exists a zero mode solution, ip ~ e~~f w , with boundary 
conditions that are given trivially by the equations of motion. In the limit fiz 0, however, 
this is equivalent to the boundary condition (|47p . Since the potential in (HT1) is equivalent to 
( 1541) up to a constant shift of the reference potential, we can conclude then that the solution 
i\) ~ e~f w is in fact a good approximation for and that, 

ml ~ (-1 - 2c)6/i 2 . (55) 
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This is clearly only valid when c < —1/2. We have checked that the expressions (IBTj) and 
(|55|) match the exact results well in this region. The full spectrum is plotted in Fig. [TJ The 
distinctive KK tower structure of the split case suggests the possibility of novel KK physics 
unlike that found in hard-wall models and may be interesting to study in other soft-wall 
bulk Higgs models as well. 



3.2 Comparison with Perturbative Expansions 

Recently, the possibility of modeling fermions by introducing non-constant bulk Dirac mass 
terms has been considered in Ref . [7J . For a single generation setup with quadratic bulk mass 
terms, the equations of motion are the same as (TT8]) . but now the mixing matrix (fl6l) can be 
written effectively as, 

-a f clk + clk^z 2 bkfi 2 z 2 \ . , 

M ~ 6 V bk ^ z2 c R k + c R k V 2 Z 2 ) ' 1 ' 

where c° L R , c l L R are constant coefficients. The effect of this non-constant bulk mass is that 
normalizable zero-modes persist (depending on the choice of the signs of c\ R ) even in the 
limit b —>■ 0. For small values of b, the spectrum may be found by treating the bulk Yukawa 
interaction as a perturbation on the 6 = solutions. 

Such an approach can be related to ours in some cases. For example, in the case of 
degenerate constant mass pieces, c° L = c° R , global unitary transformations may still be used 
to diagonalize the mass matrix when the bulk masses have the same functional form as the 
Higgs VEV. Thus, the introduction of non-constant bulk masses can be viewed as effectively 
changing the boundary conditions on the fields in such cases. 

It is interesting to note that the case considered in [7J is similar to our "split" case. In 
particular, they examine c\ = —c R and c° L = —c° R in detail. In the slightly different split 
configuration for the constant pieces of the bulk mass, c° L = ±1 — c° R , analytical solutions 
can be obtained in a similar fashion to our earlier analysis. For c° L = —1 — c° R , we find the 
lowest lying KK tower to be: 

5 c|_ 5 o i 2 2 y (57) 
4 2 46^ 2 2 L ' ^ J ' K J 













JR— : 












= o, 



(5? 



S' + ^/S + c^VS), (00) 



(60) 

where we have defined effective parameters, b 2 = {c l L ) 2 + b 2 , and = — c\ (2c° L + 1) to 
make the comparison with (!36|) - (l39|) clear. Note that there remains a lower bound on the 
mass for c° L 3> 1/2. We have checked that this solution describes the large c° L behavior for 
the case considered in Ref. [TJ, c° L = —c R . We expect that all of the basic features of the 
non-constant bulk mass model should be contained within our exact solutions. 
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Generically the split and degenerate cases allow one to find the spectrum exactly by 
solving a set of decoupled second-order equations. Even when the equations cannot be solved 
exactly, approximate methods such as those we have described above may be employed. 
Additionally, as our numerical results will verify, one can expect the behavior in these special 
cases to provide a complete qualitative picture of the full parameter space dependence. 

3.3 Couplings to Gauge Bosons 

Of significant interest in models involving extra dimensions is the coupling of fermions to 
the KK gauge bosons. When the fermions are localized at different points along the ex- 
tra dimension, they can obtain non-universal couplings to the excited gauge bosons. Such 
non-universality will generically lead to large contributions to flavor physics observables, 
providing very stringent lower bounds on the allowed KK scale [14J. 

In hard-wall models, the couplings can become universal for certain regions of the pa- 
rameter space, resulting in a GIM-like suppression of flavor changing neutral currents [3], 
thereby greatly lowering the bound on the allowed KK scale. We therefore would like to see 
if a similar effect is present in the soft-wall case. Moreover, we would like to develop our 
formalism in such a way that multiple fermion generations can be incorporated. 

The couplings of the zero mode fermions to the KK gauge bosons are found to be: 



where f A is the gauge boson profile along the extra dimension. The gauge boson profiles 
arising from a quadratic dilaton ([3]) were derived in [5] . Here we simply use the results. The 
zero-mode couplings g^^ = g5 al3 remain universal due to the orthonormality condition 
(JTTJ) and the flat zero mode gauge boson profile. This is because the dilaton factor explicitly 
cancels and plays no role. 

The degenerate single-generation case was considered in [Sj , where it was found that only 
one of the couplings, g + or g_, can become universal due to the opposite localizations of 
the fermion modes. In Ref [7], it was seen that opposite constant and non-constant bulk 
masses led to universal couplings for both g + and g_. This happens as well for the split case 
solutions. We have plotted the couplings for this case in Figs. |2] and [31 We find that the 
couplings g + and g_ become universal simultaneously whenever c^> 1/2. 

Note that the bounds from flavor physics are generically expected to be more stringent 
in soft-wall models than in models with a hard wall. This follows from the generically 
closer spacing of the KK modes in soft-wall models as compared to hard-wall models. For 
example, we can consider the contribution to ArriK arising from non-universal couplings. 
The effective 4D Lagrangian contains operators that are suppressed by the squared masses 
of the KK gauge bosons mediating the strangeness-changing transitions (AS* = 2): 




(61) 




(62) 
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Figure 2: The ratio g\j g for n = 1 (solid), n = 2 (dashed), and n = 3 (dotted) KK gauge 
modes coupled to the zero-mode fermion in the split calculated using fl6T|) . 




Figure 3: The ratio g™/g for n — 1 (solid), n = 2 (dashed), and n = 3 (dotted) KK gauge 
modes coupled to the zero-mode fermion in the split CclS6 ; clS calculated using (|6ip . 
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where the sum is over the gauge boson KK modes with KK masses M n , and g± = 
R 9± V^ R with R generic unitary matrices [14J . Thus, in the presence of non- 
degenerate couplings to the bulk KK gauge bosons, bounds from flavor experiments may 
be interpreted lower bound on the KK scale. 

The key point is that the total amount of suppression in (1621) depends upon the spacing 
of the KK tower. In a hard- wall model, for example, m? n ~ n 2 M^- K , where Mkk is the 
KK mass scale. This compares with the soft-wall scenario where it would seem to imply a 
problem, because the squared mass trajectories grow generically as rv? n ~ nM\ K (indeed, 
this spacing was the original motivation for studying the soft-wall While the sum of 

1/n diverges as n — * oo, we should of course truncate the sum at some high energy cutoff. 
Nevertheless, the naive implication is that the constraints on soft-wall models should be 
considerably tighter. 

However, this argument ignores the fact that the gauge bosons become increasingly IR 
localized with increasing mode number n. Thus, any off-diagonal terms in the gauge coupling 
matrices are further suppressed for large n. By performing a numerical fit using the first 
several dozen gauge boson modes and our split case solutions, we find that the couplings fall 
off as n~ 0A to a very good approximation in the region where the couplings are independent 
of localization. This implies that the terms in the sum (162|) grow as n~ 1A . All other things 
being equal, this implies that the constraints from flavor physics are roughly a factor of two 
more stringent in this model than in hard-wall models. 

While this presents no great problem for the model with a quadratic dilaton, for a 
generic power law behavior in the dilaton, $ ~ z a , the spectrum of gauge bosons grows 
as w? n ~ n 2 ~ 2 / a Mx K [5J. This means that for less steep potentials, even tiny amounts 
of non- degeneracy among the bulk couplings has potentially severe implications for flavor 
physics. 

4 Numerical Solution 
4.1 Routine 

The analytical solutions that we have presented are of limited use, and instead we would like 
to solve the full fermion mass problem including flavor. Our goal is to solve the eigenvalue 
problem ffT5]) with mixed boundary conditions. The "initial conditions" ([5]) specify half of 
the boundary values at the UV brane or, equivalently, half of the integration constants for 
the system. The remaining constants of integration are fixed by the normalization conditions 
( TTTj) . which can only be satisfied if the eigenvalue, m", has been chosen correctly. 

We convert the problem to an initial value one by extending the shooting method to 
linear boundary value problems of arbitrary order [16J. The solutions to ffTol) may be written 
as: 

r (n) {z) = U{m a. ^ Zo yJfMn) {zo)j (63) 

where the propagator U(m";z, z$) is a linear operator and the f %a(jL \z) are ANp x ANp 
matrix-valued functions for Np fermion generations. 
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log 10 (mo/VzT jJ.) 
lr 

split 




Figure 4: Lowest lying masses for the "degenerate" (c = cl = cr) and "split" (c = cl = 
— 1 — cr) cases, where the solid lines are determined in Ref. [5] (degenerate) and from (148!) 
(split). The dots represent values obtained using the numerical method of Section I4TT1 

The matrix elements of U may be found by integrating a set of 4Np linearly independent 
basis vectors that span the space of initial values, f %a ^(zo), and inverting (1631) . The initial 
values that lead to normalizable solutions correspond to eigenvectors of U(m"; z, z ) with 
vanishing eigenvalues in the limit z — > oo. There are generally 2Np such eigenvectors. 
Numerically, we can estimate the values of these vectors by considering the eigenvectors 
of U(m%; zi, Zo), where our cutoff satisfies z\ ^> fi^ 1 . In practice, results are much more 
reliable if one starts the forward integration from some intermediate range z* ~ fi^ 1 and 
then integrates the normalizable modes back to Zq. Variations on this theme can be explored. 

We scan over m®, at each point integrating the system from a set of initial values so as to 
reconstruct the normalizable solutions. If there exists a linear combination of the solutions 
that matches the boundary conditions (j5J), then is a solution to the system. To determine 
when this occurs, we define a merit function as the absolute value of the determinant of a 
matrix and search for a minimum. The matrix we use has columns formed by projecting 
out of the normalizable initial value vectors those components that are not restricted by the 
boundary conditions. These projected column vectors must be linearly dependent in order 
to satisfy the boundary conditions of the problem. 

When the hierarchy between fj, and Zq 1 is very large, increasingly high precision is nec- 
essary to achieve reliable results. Iterative methods may be better suited to the problem in 
such cases. Our primary goal is to highlight the differences between fermions in soft-wall and 
hard-wall scenarios, and the speed and simplicity of this technique are its chief advantages. 
For this reason, we have limited our attention to a modest hierarchy. 
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Figure 5: Contours of log 10 (m / \/6a0 for the lowest lying masses in our soft-wall setup with 
Vbfjizo = 10~ 3 . 

4.2 Results 

4.2.1 Single Generation 

We first present results for a single generation of fermions, as this case illustrates the essential 
features of the fermion mass behavior in the soft- wall, and allows us to compare our numerical 
results with the analytical cases in the appropriate limits as well as to a typical hard-wall 
setup. Assume the following values of the parameters: 

H = 1 TeV; fiz = 10~ 3 ; 6=1. (64) 

In Figure HI we compare the numerical results to the analytical results from Section 13.11 
where it can be seen that the two methods agree very well. Next in Figure [5] we plot the 
fermion mass contours to show the full dependence on the parameters cl and cr . The shape 
of the plot is easily understood from the analytical results. The numerical solution smoothly 
interpolates between the solutions along the lines cl = cr and cx, = ±1 — c R . Because a 
similar analysis can be repeated for other Higgs VEVs, this provides a natural way to begin 
studying the qualitative aspects of other models in AdS as well. 

We can compare the soft-wall behavior with a typical hard-wall setup. In Figure [6] we 
provide the corresponding contour plot for a hard- wall model in which the SM fermion masses 
are simply proportional to the values of the wavefunctions at z = l//i. The most striking 
difference between the plots occurs in the lower right-hand corner. This is the region where 
c L > 1/2 and c R < -1/2. 
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Figure 6: Contours of log 10 (m /\/&/i) for the lowest lying masses in a typical hard- wall setup 
with Vbfiz = 1CT 3 . 

The hard-wall case is characterized by a steep dependence on the bulk mass in this region, 
where the wavefunctions are proportional to ffj_ ~ z~ CL and ~ z CR . For zq <C the 
normalization constants become vanishingly small: 

Afi^o" 1/2±C ^ (65) 

Thus, the values of the functions in the IR at z = are additionally suppressed. This is 
the well-known mechanism for generating SM mass hierarchies in Randall-Sundrum scenarios 
with bulk fields [21 Ej. For the soft- wall case, however, we can see the lower bound on the 
mass in this region, 

m ~ OuzoV, (66) 

as indicated by the approximate expression fl5Tj) . This can be understood be noting that 
the normalization ffl7|) involves the sum of two types of fermion contributions which are 
generically not simultaneously suppressed. 

4.2.2 Three Generations 

Next we aim to provide concrete numerical examples involving three generations of fermions 
that fully take into account the 5D flavor mixing to show that the attractive features of the 
soft-wall are maintained. For multiple generations, there are three matrices that parame- 
terize the fermions: two bulk mass matrices and M#, and the bulk Yukawa matrix, A5. 
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We take the action (j4]) to be written in an arbitrary basis, for example, the CKM basis. 
Absent some symmetry, there is no reason to expect any structure relating the entries of the 
various bulk parameter matrices. We generically expect that the entries of each matrix are 
all of order unity (in units of the AdS curvature scale, k), and that the various matrices are 
misaligned. There is of course some basis in which both Ml and Mr are diagonal. Thus, by 
"misaligned," we mean that this basis is distinct from the one in which the Yukawa matrix 
is diagonal. Indeed, the typical approach is to work in this basis, treating the Yukawa inter- 
actions as perturbations. Such an approach has been used in both hard- wall [2j [3J [T71 [8j [18] 
and even very recently in soft- wall setups [Tj. 

In Ref.jTj, it was found that one needed to include the first several (~ 10) KK modes in 
order to achieve reliable results in such a perturbative expansion when including only a single 
generation. At such a point, the analysis is essentially a numerical exercise. In our view, 
it is advantageous to include the entire KK tower in the numerical formulation wherever 
possible. In other words it may make the most sense to simply solve the equations of motion 
f ll5p . which guarantee the orthogonality of the eigenfunctions due to the hermiticity of the 
mixing matrix. 

We expect that all other interactions may be treated reliably as perturbations. This 
is because the Higgs grows unbounded in the IR where it is the dominant contribution to 
the fermion equations of motion. Other observables may thus be calculated using the usual 
wavefunction overlap approximation. As an application, we will calculate the couplings to 
excited gauge bosons for examples involving three generations. 

We do not attempt to set precise bounds on soft-wall models here, as doing so goes 
significantly beyond the scope of this work. Electroweak and flavor constraints have been 
discussed in the context of soft-wall models in Refs. [3 [7]. Detailed analyses in various 
hard- wall scenarios can be found in [HI [2J El EH E333 EES] and references therein. 

However, we will require that the eigenvalues of the bulk mass matrices satisfy m l L > k/2 
and m l R < — m l L in order to get nearly degenerate gauge couplings. Because of the lower 
bound on the fermion masses at mo ~ (fizo)[i in this region, it is clear that the hierarchy 
considered above, fizo = 10~ 3 will be inadequate for generating MeV scale masses when 
\i = 1 TeV, and will only be possible for fiz < 10~ 6 . Thus we again assume a quadratic 
Higgs VEV, h(z) = rik 3 / 2 fj?z 2 , and the following for our input parameters: 



Dealing with much larger hierarchies presents significant numerical challenges. However, the 
qualitative results of such an analysis should not be substantially different from the results 
presented here. 

First, we present an example resembling down- type quarks (or charged leptons). For 
simplicity, we take the entries of Mi to be nearly degenerate, but we allow for large non- 
degeneracy in the matrix M# as well as in the Yukawa matrix. Specifically, we consider, 



1 TeV; 



jiz 



10 



6 



(67) 
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Figure 7: The down- type fermion bulk profiles F"(z) (in units of ^fjX) for the first generation 
(solid), second generation (dashed) and third generation (dotted) showing the overlap with 
the Higgs VEV h(z) (in units of jj? / \fk with rj = 1). 



/ 0.422 0.175 -0.678 \ 
-^=A 5 = -0.007 0.928 0.348 . (68) 
\ 0.295 -0.327 0.637 / 

We find a spectrum of masses resembling the down-type quarks (or charged leptons): 

m£ = 0.57 MeV, 96.08 MeV, 1.310 GeV. (69) 

The fermion mass hierarchy is clearly obtained, but due to the complexity of the numerical 
procedure we do not match the SM masses exactly, and postpone a more detailed analysis 

for future work. The fermion bulk profiles, F®(z) = J {f^-^y Tl-^ + {fR-^ Tr-^ are 
plotted in Figure [71 The fermion profile overlap with the Higgs VEV, h(z) leads to the 
fermion mass hierarchy. The corresponding bulk profiles, F"(z) are not plotted because the 
profile differences between the flavors are not as pronounced. This is due to our choice of 
UV boundary conditions and bulk masses (1681) . 

From expression (IBT]) . we can calculate the coupling of the zero mode fermions to the 
KK gauge bosons (i.e. gluons). The result is a matrix whose off-diagonal entries contribute 
to flavor violation. We obtain the following results for the first two KK gauge coupling 
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Figure 8: The up- type fermion bulk profiles F"(z) (in units of yfjx) for the first generation 
(solid), second generation (dashed) and third generation (dotted) showing the overlap with 
the Higgs VEV h(z) (in units of jj? / \fk with rj = 1). 



matrices, normalized to the coupling to the massless gauge boson: 
(i) / 0.186 10" 4 1(T 4 \ ( 2 ) 

_ I in-4 o lg7 2 x 1Q -4 I '/+ 

2 x 10~ 4 0.185 




0.188 ^0 ^0 \ ( 2 ) 
38 10~ 4 , 9 — 
10~ 4 0.184 / 9 




0.188 KT 4 I , — = I « 0.139 10" 4 I . (70) 



This behavior is maintained for higher modes as well. For this choice of parameters, the very 
nearly degenerate couplings imply that \i of order a few TeV will be consistent with flavor 
constraints [I21EIE1E]- Note that we have assumed no contributions to CP violation. Thus 
the soft-wall model can accommodate the fermion mass hierarchy with large bulk mixing 
and small flavor violation. 

The up-type quarks are only moderately more sensitive to the presence of the top quark 
when large bulk mixing is allowed. For the choices (|67|) we obtain 

n , ( 0.749 -0.005 0.017 \ / -0.940 -0.285 -0.200 

Ml / „ I Mr 



0.005 0.785 0.066 , ~T = \ -0-285 -1.103 -0.338 
0.017 0.066 0.516 / \ -0.200 -0.338 -0.657 



0.700 -0.352 -0.193 
-0.079 0.826 -0.065 
-0.098 -0.321 1.430 



-^=A 5 = | -0.079 0.826 -0.065 | , (71) 
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which gives rise to the following mass spectrum: 



m° = 2.10 MeV, 129.1 MeV, 151.5 GeV. (72) 

Again we see that the correct fermion mass hierarchy can be obtained. The fermion bulk 
profiles, F®(z) are plotted in Figured The fermion profile overlap with the Higgs VEV, h(z) 
leads to the fermion mass hierarchy. Similarly to the down-type fermions, the corresponding 
up- type bulk profiles F+(z) are not plotted because the profile differences are negligible due 
to the choice of UV boundary conditions and bulk masses (1711) . The gauge couplings are 
nearly universal among the first two generations: 

0.140 10~ 3 10~ 3 \ 
10- 3 0.139 10- 3 ; 
10- 3 10- 3 -0.05 / 

0.140 10~ 6 3 x 10" 4 \ 
lO" 6 0.137 2 x lO" 3 . 
3 x 10~ 4 2 x 10~ 3 -0.140 / 

(73) 

Constraints from top quark physics are significantly weaker, so this is not expected to affect 
the bound on \i. 



9 



9 



(i) 
+ 

9 

(i) 



9 




2 x lO" 3 2 x 10~ 3 

0.185 10~ 6 

10- 6 -0.05 

2 x 10^ 6 4 x 10" 4 

0.183 10~ 3 

10~ 3 -0.17 



(2) 

9+ 



9 



9 

(2) 

9 



5 Summary 

We have presented a variety of tools useful for studying fermion physics in soft-wall back- 
grounds, focusing heavily on the treatment of fermion masses. The equations of motion are 
non-trivial to solve and generically require numerical techniques. However, we have doc- 
umented several special cases for which it is possible to decouple the equations of motion. 
These cases serve as useful examples for qualitatively understanding the full parameter space 
behavior, as they illuminate independent "axes" of the parameter space along which fermion 
behavior can be understood in detail. The utility of our approach is due not only to the fact 
that it effectively reduces the problem to solving a set of one-dimensional Schrodinger-like 
equations, for which many theoretical and numerical tools have been created, but also to 
the fact that it applies to any soft-wall model in AdS space. This opens up the possibility of 
analyzing fermions in a wide variety of Higgs models. For example, it should be possible to 
analyze fermion physics in unHiggs scenarios, such as that considered in Ref. [9], or to study 
other power-law Higgs behavior, as has been examined in the degenerate case in Refs. [51 E] ■ 
Furthermore, we have outlined methods for calculating fermion masses and wavefunc- 
tions in an arbitrary background with multiple flavors and arbitrary bulk parameters. The 
formalism maintains the orthogonality of the KK tower, making it particularly useful for 
studying the experimental consequences of soft-wall models with additional bulk fields. For 
example, we showed explicitly how to calculate the fermionic couplings to KK gauge bosons. 
The off-diagonal entries in the coupling matrix are directly related to the amplitudes of flavor 
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changing neutral current processes. Moreover, we argued that the experimental constraints 
on new sources of flavor violation are generically tighter in soft-wall models than in hard-wall 
models due to the smaller spacing of the KK resonances. While the tightening is not too 
constraining in a model with a quadratic dilaton, for an IR cut-off growing much less quickly 
than z 2 , the constraints can become severe for model building. 

We described a very simple numerical technique for calculating fermion spectra, which 
we used to show the full parameter space dependence of SM fermion masses on the bulk 
mass parameters. The technique maintains the attractive features of the formalism, such 
as an orthogonal KK tower and extends naturally to incorporate several generations of 
fermions. Thus, the solutions allow for the straightforward calculation and interpretation of 
new physics observables. 

Both the analytical and numerical results suggest the potential for rich collider physics 
that is substantially different from that obtained in hard-wall models. While our results were 
for a particular choice of a bulk Higgs VEV, they demonstrate that a soft-wall background 
can lead to a distinctive phenomenology. For example, our results indicate the presence of 
a lower bound on fermion masses in a large area of the parameter space, suggestive of a 
seesaw-like mechanism. This could easily be implemented to explain neutrino masses in a 
scenario where the hierarchy between k and \x is of order the GUT scale. By introducing three 
additional right handed neutrino fields, a fairly random difference between the bulk masses 
of the right-handed neutrinos and charged leptons could naturally lead to light neutrino 
masses. This would be an interesting variation on the ideas that are well-known in the 
hard- wall picture (cf. QH [20] ) . 

Even with these differences, the essential and attractive features of the hard-wall can be 
retained in our model. We presented results for three generations of fermions with anarchic 
5D parameters that reveal standard model-like particle masses and GIM-like suppression of 
KK gauge boson mediated flavor changing neutral currents. We argued that this implied a 
fairly modest bound on the KK scale. A more general analysis of flavor physics bounds will 
lead to stringent constraints in the soft-wall model. Having developed the tools needed to 
examine electroweak and flavor physics in full detail, a more detailed study can be now be 
undertaken. 
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